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Abstract. To any complex Hadamard matrix we associate a quantum permu- 
tation group. The correspondence is not one-to-one, but the quantum group 
encapsulates a number of subtle properties of the matrix. We investigate various 
aspects of the construction: compatibility to product operations, characterization 
of matrices which give usual groups, explicit computations for small matrices. 



Introduction 

A complex Hadamard matrix is a matrix h G M n (C), having the following prop- 
erty: entries are on the unit circle, and rows are mutually orthogonal. Equivalently, 
n~ l l 2 h is a unitary matrix with all entries having the same absolute value. 

The basic example is the Fourier n x n matrix, given by F^s = w y , where w = 
e 2 W n . There are many other examples, but there is no other known family of 
complex Hadamard matrices that exist for every n. 

These matrices are related to several questions: 

(1) The story begins with work of Sylvester, who studied real Hadamard matrices 
[5]. These have ±1 entries, and can be represented as black and white 
pavements of a n x n square. Such a matrix must have order n = 2 or 
n = 4k, and the main problem here is whether there are such matrices for 
any k. At the time of writing, the verification goes up to k = 166. 

(2) One can try to replace the ±1 entries by numbers of the form ±1, ±i, or by 
n-th roots of unity (the Fourier matrix appears), or by roots of unity of ar- 
bitrary order. These are natural generalizations of real Hadamard matrices, 
first investigated by Butson in the sixties |Buj . 

(3) The matrices with arbitrary complex entries appeared in the eighties. Popa 
discovered that such a matrix gives rise to the operator algebra condition 
A _L hAh* , where A C M n (C) is the algebra of diagonal matrices |Poj . Then 
work of Jones led to the conclusion that associated to h are the following 
objects: a statistical mechanical model, a knot invariant, a subfactor, and 
a planar algebra. The computation of algebraic invariants for these objects 
appears to be remarkably subtle and difficult. See Q]. 

(4) In the meantime, motivated by a problem of Enffo, Bjorck discovered that 
circulant Hadamard matrices correspond to cyclic n-roots, and started to 
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construct many examples [Bj]. In a key paper, Haagerup proved that for 
n = 5 the Fourier matrix is the only complex Hadamard matrix, up to per- 
mutations and multiplication by diagonal unitaries jHj- Self-adjoint complex 
Hadamard matrices of order 6 have been recently classified BN . 

(5) Another important result is that of Petrescu, who discovered a one-parameter 
family at n = 7, providing a counterexample to a conjecture of Popa regard- 
ing the finiteness of the number of complex Hadamard matrices of prime 
dimension |Pej . The notion of deformation is further investigated in |N] . 

(6) Several applications of complex Hadamard matrices were discovered in the 
late nineties. These appear in a remarkable number of areas of mathemat- 
ics and physics, and the various existence and classification problems are 
fast evolving. See |TZj for a catalogue of most known complex Hadamard 
matrices and a summary of their applications. 



The purpose of this paper is to present a definition for what might be called 
symmetry of a Hadamard matrix. This is done by associating to h a certain Hopf 
algebra A. This algebra is of a very special type: it is a quotient of Wang's quantum 
permutation algebra [Wa2j . In other words, we have the heuristic formula A = C(G), 
where G is a quantum group which permutes the set {1, . . . , n}. 

As an example, the Fourier matrix gives G = 7L n . In the general case G doesn't 
really exist as a concrete object, but some partial understanding of the quantum 
permutation phenomenon is available, via methods from finite, compact and discrete 
groups, subfactors, planar algebras, low-dimensional topology, statistical mechanical 
models, classical and free probability, geometry, random matrices etc. See BB , 
BBC], [BClj . |BC2j and the references there. The hope is that the series of papers on 
the subject will grow quickly, and evolve towards the Hadamard matrix problematics. 
A precise comment in this sense is made in the last section. 

Finally, let us mention that the fact that Hadamard matrices produce quantum 
groups is known since |Balj . |Wa2j . The relation between coinvolutive compact 
quantum groups, abstract statistical mechanical models and associated commuting 
squares, subfactors and standard A-lattices is worked out in |Ba2j . Another result is 
the one in |Ba3j . where a Tannaka-Galois type correspondence between spin planar 
algebras and quantum permutation groups is obtained. In both cases the situations 
discussed are more general than those involving Hadamard matrices, where some 
simplifications are expected to appear. However, the whole subject is quite technical, 
and it is beyond our purposes to discuss the global picture. This paper should be 
regarded as an introduction to the subject. 

The paper is organized as follows. 1 is a quick introduction to quantum permuta- 
tion groups. 2 contains a few basic facts regarding complex Hadamard matrices, the 
construction of the correspondence, and the computation for the Fourier matrix. In 
3-6 we investigate several aspects of the correspondence: compatibility between ten- 
sor products of quantum groups and of Hadamard matrices, the relation with magic 
squares, characterization of matrices which give usual groups, explicit computations 
for small matrices, and some comments about subfactors and deformation. 
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1. Quantum permutation groups 

This section is an introduction to A(S n ), the algebra of free coordinates on the 
symmetric group S n . This algebra was discovered by Wang in |Wa2j . 

The idea of noncommuting coordinates goes back to Heisenberg, the specific idea 
of using algebras of free coordinates on algebraic groups should be attributed to 
Brown, and a detailed study of these algebras, from a K-theoretic point of view, is 
due to McClanahan. Brown's algebras are in fact too big, for instance they have 
no antipode, and the continuation of the story makes use of Woronowicz's work 
on the axiomatization of compact quantum groups. The first free quantum groups, 
corresponding to U n and O n , appeared in Wang's thesis. The specific question about 
free analogues of S n was asked by Connes. See |Brj . [Mj . |Walj . [Wa2 1- 

We can see from this brief presentation that the algebra A(S n ) we are interested in 
comes somehow straight from Heisenberg, after a certain period of time. In fact, as 
we will see now, the background required in order to define A(S n ) basically reduces 
to the definition of C*-algebras, and to some early work on the subject: 

Let A be a C*-algebra. That is, we have a complex algebra with a norm and an 
involution, such that Cauchy sequences converge, and ||oo*|| = ||o|| 2 - 

The basic example is B(H), the algebra of bounded operators on a Hilbert space 
H. In fact, any C*-algebra appears as subalgebra of some B(H). 

The key example is C(X), the algebra of continuous functions on a compact space 
X. This algebra is commutative, and any commutative C*-algebra is of this form. 

The above two statements are results of Gelfand-Naimark-Segal and Gelfand. 
The proofs make use of standard results in commutative algebra, complex analysis, 
functional analysis, and measure theory. The relation between the two statements 
is quite subtle, and comes from the spectral theorem for self-adjoint operators. The 
whole material can be found in any book on operator algebras. 

Wang's idea makes a fundamental use of the notion of projection: 

Definition 1.1. Let A be a C* -algebra. 

(1) A projection is an element p € A satisfying p 2 = p = p* . 

(2) Two projections p,q E A are called orthogonal when pq = 0. 

(3) A partition of unity is a set of orthogonal projections, which sum up to 1. 

A projection in B(H) is an orthogonal projection P(K), where K C H is a closed 
subspace. Orthogonality of projections corresponds to orthogonality of subspaces, 
and partitions of unity correspond to decompositions of H. 

A projection in C(X) is a characteristic function x(X)i w here Y C X is an 
open and closed subset. Orthogonality of projections corresponds to disjointness of 
subsets, and partitions of unity correspond to partitions of X. 
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Definition 1.2. A magic unitary is a square matrix u £ M n {A), all whose rows 
and columns are partitions of unity in A. 

The terminology comes from a vague similarity with magic squares, to be investi- 
gated later on. For the moment we are rather interested in the continuing the above 
classical/quantum analogy, for projections and partitions of unity: 

A magic unitary over B{H) is of the form P(Kij), with K magic decomposition 
of H, in the sense that all rows and columns of K are decompositions of H. The 
basic examples here are of the form = C£j,-, where £ is a magic basis of H, in 
the sense that all rows and columns of £ are bases of H. 

A magic unitary over C(X) is of the form x(Yij), with Y magic partition of X, in 
the sense that all rows and columns of Y are partitions of X. The key example here 
comes from a finite group G acting on a finite set X: the characteristic functions 
Xij = {a £ G \ a(j) = i} form a magic unitary over C(G). 

In the particular case of the symmetric group S n acting on {1, . . . , n}, we have 
the following presentation result, which follows from the Gelfand theorem: 

Theorem 1.1. C(S n ) is the universal commutative C* -algebra generated by n 2 ele- 
ments Xij > with relations making (xij) a fno,gic unitary matrix. Moreover, the maps 



S(Xij) = Xji 
are the comultiplication, counit and antipode of C{S n ). 

In other words, when regarding S n as an algebraic group, the relations satisfied 
by the n 2 coordinates are those expressing magic unitarity. Indeed, the characterstic 
functions Xij are nothing but the n 2 coordinates on the group S n C O n , where the 
embedding is the one given by permutation matrices. 

See the preliminary section in [B"B~C] for a proof of the above result, and the 
introduction of BC2 for more comments on the subject. 

We are interested in the algebra of free coordinates on S n . This is obtained by 
removing commutativity in the above presentation result: 

Definition 1.3. A(S n ) is the universal C* -algebra generated by n 2 elements x^, 
with relations making (xij) a magic unitary matrix. The maps 



are called the comultiplication, counit and antipode of A{S n ). 

This algebra, discovered by in |Wa2j , fits into the quantum group formalism devel- 
oped by Woronowicz in |Wolj . |Wo2j . In fact, the quantum group G n defined by the 
formula A(S n ) = C(G n) is a free analogue of the symmetric group S n . This quan- 
tum group doesn't exist of course: the idea is just that various properties of A(S n ) 
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can be expressed in terms of it. As an example, the canonical map A(S n ) — ► C(S n ) 
should be thought of as coming from an embedding S n CZ Gn- 

Proposition 1.1. For n = 1,2,3 we have A(S n ) = C(S n ). 

This follows from the fact that for n < 3, the entries of a n x n magic unitary 
matrix have to commute. For instance at n = 2 the matrix must be of the form 

\l-p p J 

where p is a projection, and entries of this matrix obviously commute. 

The result is no longer true for n = 4, where more complicated examples of magic 
unitary matrices are available, for instance via diagonal concatenation: 



/ 


p 


1 — p 








\ 


1 


-p 


p 


















q 


1 - 


q 


V 








l-q 


Q 


J 



In fact, A(S n ) with n > 4 is not commutative, and infinite dimensional. 
Consider now an arbitrary magic unitary matrix u £ M n (A). We say that it is a 
corepresentation of A if its coefficients generate A, and if the formulae 

A(Uy) = ^ u ik ® u k j 

e(uij) = Sij 
S{u,ij) = Uji 

define morphisms of C*-algebras. These morphisms have by definition values in the 
algebras A (g> A, C and A op , and they are uniquely determined if they exist. In case 
they exist, these morphisms make A into a Hopf algebra in the sense of E Wol . 
This notion provides an axiomatization for quotients of A(S n ): 

Definition 1.4. A quantum permutation algebra is a C* -algebra A, given with a 
magic unitary corepresentation u 6 M n (A). 

We have the following purely combinatorial approach to these algebras: 
Recall first that the unitary representations of a discrete group L are in correspon- 
dence with representations of the group algebra C*(L). The same is known to hold 
for discrete quantum groups, so we have a correspondence between representations 

L n -> U{H) 
A(S n ) - B{H) 

where L n is the discrete quantum group whose group algebra is A(S n ). We can 
therefore consider the quantum permutation algebra A = C* (Im(L n )). 
All this is quite heuristic, but the construction of A is definitely possible: 

Definition 1.5. Associated to a representation ir : A(S n ) — > B(H) is the smallest 
quantum permutation algebra A n realizing a factorization of n. 
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In this definition we assume of course that the morphism making tt factorize is 
the canonical one, given by Xij — > Uij. The construction of A n is standard, for 
instance by taking the quotient of A(S n ) by an appropriate Hopf algebra ideal. The 
uniqueness up to isomorphism is also clear. See Bal . 

We have the following examples: 

(1) For the counit representation tt : A(S n ) ->Cwe get A^ = C. 

(2) For a faithful representation tt : A(S n ) C B(H) we get A n — A(S n ). 

(3) If A is a quantum permutation algebra, A C B(H), and tt : A(S n ) — ► A is a 
surjective representation of A(S n ), then A n = A. 

These examples are all trivial. Note however that the third one has the following 
interesting consequence: 

Proposition 1.2. Any quantum permutation algebra is of the form A n . 

Proof. Let A be a quantum permutation algebra. We can compose the canonical 
quotient map A(S n ) — > A with an embedding A C B(H), say given by the GNS 
theorem, and we get a representation tt as in the statement. □ 

2. Hadamard matrices 

Let h G M n (C) be an Hadamard matrix. This means that all entries of h have 
modulus 1, and that rows of h are mutually orthogonal. In other words, we have 

( kl \ 

h= h2 

\k) 

where the vectors hi are formed by complex numbers of modulus 1, and are orthog- 
onal with respect to the usual scalar product of C n , given by: 

< x, y >= xiyi + . . . + x n y n 

It follows from definitions that the columns of h are orthogonal as well. 
We have the following characterization of such matrices: 

Proposition 2.1. For a matrix h £ M n (C), the following are equivalent: 

(1) h is an Hadamard matrix. 

(2) n~ l / 2 h is a unitary matrix, all whose entries have same modulus. 

We should mention that in the operator algebra literature, it is rather n _1 / 2 /i 
that is called Hadamard. The other warning is that in the combinatorics literature, 
the Hadamard matrices are just the real ones, meaning those having ±1 entries. We 
hope that the hybrid terminology used in this paper won't cause any trouble. 

We use capital letters to denote explicit Hadamard matrices. The first example, 
which is in fact the only basic example, is the Fourier matrix: 

Definition 2.1. The Fourier matrix of order n is given by Fij = w 13 , where w = 

p 2iri/n 
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We have a natural equivalence relation for Hadamard matrices, given by h ~ k if 
one can pass from h to k by permutations of rows and columns, and by multiplica- 
tions of rows and columns by complex numbers of modulus 1. See |TZ| . 

For n= 1,2,3,5 any Hadamard matrix is equivalent to the Fourier one, see |H| . 
At n = 4 we have the following example, depending on q on the unit circle: 



M q = 



fl 1 1 


1\ 


1 q -1 


-Q 


1 -1 1 


-1 


\1 ~Q "I 


9/ 



These are, up to equivalence, all 4 x 4 Hadamard matrices. As an example, the 
Fourier matrix corresponds to the value q = ±i. See Haagerup [Hj . 
At n = 6 we have several examples, as for instance: 



H = 



fi 


1 


1 


1 


1 


1\ 


1 


i 


1 


-1 


-1 


1 


1 


1 


i 


1 


-1 


-1 


1 


-1 


1 


i 


1 


-1 


1 


-1 


-1 


1 


i 


1 


V 1 


1 


-1 


-1 


1 


4 / 



This matrix appears in the paper of Haagerup [H]. Another example appears in 
the paper of Tao [T]. This is given by the following formula, with w = e 27 ™/ 3 : 



T = 



fl 


1 


1 


1 


1 


1\ 


1 


1 


«' 




w 2 


w 2 


1 


«' 


1 


w 2 


UJ 2 


u; 


1 


«' 


?/r 


1 


w 


u; 2 


1 


w 2 






1 




V 


w 2 


?/; 


u; 2 


?(' 





Observe that the Fourier, Haagerup and Tao matrices are based on certain roots 
of unity. We have here the following notion: 

Definition 2.2. The Butson class Hi(n) consists of Hadamard matrices h € M n (C) 
having the property h\j = 1 for any i,j. 

In this definition I is a positive integer. It is convenient to use as well the notation 
Hoo{n), for the class of all n x n Hadamard matrices. Here are a few examples: 

(1) The Fourier matrix is in H n {n). 

(2) The Haagerup matrix is in H^(6). 

(3) The Tao matrix is in #3(6). 

(4) The real Hadamard matrices are in H^in). 

(5) The matrix M(q) is in i?;(4), where I is the order of q 2 . 

We call n size, and I level. In lack of some better idea, the complexity of an 
Hadamard matrix will be measured by its size and level, the size coming first. 
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At infinite level the self-adjoint matrices of order 6 are classified in |BN| . The 
complete list of Hadamard matrices at n = 6 is not known. 
We end this discussion with the following well-known fact: 

Proposition 2.2. Each Hadamard matrix is equivalent to a Hadamard matrix which 
is dephased, meaning that the first row and column and the diagonal consist of 1 's. 

Consider now an arbitrary Hadamard matrix h S M n (C). The rows of h, denoted 
as usual h%, . . . ,h n , can be regarded as elements of the algebra C n . 

Since each hi is formed by complex numbers of modulus 1, this element is invert- 
ible. We can therefore consider the following matrix of elements of C n : 

p..- h 

KtJ ~ hi 

The scalar products on rows of £ are computed as follows: 

<iij,iik> = < hj/hi,hk/hi > 
= n < hj,hk > 
= n ■ 5 jk 

In other words, each row of £ is an orthogonal basis of C n . A similar computation 
works for columns, so £ is a magic basis of C n . Thus we can apply the procedure 
from previous section, and we get a magic unitary matrix, a representation, and a 
quantum permutation algebra: 

Definition 2.3. Let h S M n (C) be an Hadamard matrix. 

(1) £(h) is the magic basis given by £ij = hj/hi. 

(2) P(h) is the magic unitary given by Pij = P(£,ij)- 

(3) Tih is the representation given by -Kh{uij) = Pij. 

(4) A(h) is the quantum permutation algebra associated to tt. 

In other words, associated to h are the rank one projections P(hj/hi), which form 
a magic unitary over M n (C). We consider the corresponding representation 

tt : A(S n ) -> M n (C) 

we say that this comes from a representation L n — > U(n) of the dual of the n-th 
quantum permutation group, and we consider the algebra A = C*{Im{L n )). 

It is routine to check that equivalent Hadamard matrices hi, /12 give the same al- 
gebra, A{h\) = A(h 2 ). This is because the standard equivalence operations, namely 
permutation of rows and columns, and multiplication or rows and columns by scalars 
of modulus 1, give conjugate magic unitaries, hence conjugate representations of 
A(S n ). As for the converse, this is expected not to hold, see the conclusion. 

As a first example, consider the Fourier matrix: 

Theorem 2.1. For the Fourier matrix F of order n we have A(F) = C(Z n ). 
Proof. We have the following computation, where p = (w, w 2 , . . . , w n ): 

Fij tr'-> => Fi = p l 
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=> P ij =P((S- i ) 

Consider the cycle c(i) = i — 1. We regard Z n = {0, 1, . . . , n — 1} as a subgroup 
of S n , via the embedding k — > c fc . Consider the following diagram: 

A(5 n ) -> M n (C) 



I T 

C(S n ) — C(Z n ) 

We define the connecting maps in the following way: 

(1) The map A(S n ) — > M n (C) is the representation 7r. 

(2) The map A(S' ra ) — > C(S n ) is the canonical one, given by — > Xij- 

(3) The map C(S n ) — > C(Z n ) is the transpose of Z n C Sn- 

(4) The map C(Z n ) -> M n (C) is given by 5 fc -» P(p k ). 
We compute the image of Xij by the third map: 

Im (Xii) = Im (x{cr | a(j) = i}) 
= X{k | c k (j) = i} 
= X{k | j -k = i} 

Thus at level of generators, we have the following diagram: 

Xij - P^- 1 ) 



I T 

Xij > fij—i 

This diagram commutes, so the above diagram of algebras commutes as well. 
Thus we have a factorization of it through the quantum permutation algebra C(Z n ). 
Moreover, this algebra is the minimal one making tt factorize, for instance because 
it is isomorphic to the image of tt. This gives the result. □ 



3. Tensor products 

Given two Hadamard matrices h € M n (C) and k € M m (C), we can consider their 
tensor product. This is an element as follows: 

h®k G M n (C) ® M m (C) 

In order to view h ® k as a usual matrix, we have to identify the algebra on the 
right with the matrix algebra M nm (C). This isomorphism is not canonical, so we 
proceed as follows: given h, k, we define h<g> k £ M nm (C) by the formula 

(h (8) k) ia j b = hijk ab 
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where double indices ia,jb £ {1, . . . , nm} come from indices i,j £ {1, . . . , n} and 
a,b £ {1, . . . m}, via some fixed decomposition of the index set. Most convenient 
here is to use the lexicographic decomposition: {1,2,..., nm} = {11, 12, ... , nm}. 
It is well-known that the above two definitions of h ® k coincide. 

Proposition 3.1. h®k is an Hadamard matrix. 

Proof. We can use for instance Proposition 2.1: since h, k are Hadamard matrices, 
n _1 / 2 /i, m~^l 2 k are unitaries. Now the tensor product of two unitaries being a 
unitary, we get that {nm)~ l l 2 h ® k is unitary, so h ® k is Hadamard. □ 

It is known from BB that various operations at level of combinatorial objects 
lead to similar operations at level of quantum permutation algebras. In this section 
we prove such a result for tensor products of Hadamard matrices: 

Theorem 3.1. We have A{h ® k) = A(h) ® A(k). 

Proof. It is convenient here to use tensor product notations, at the same time with 
bases and indices. We use the lexicographic identification of Hilbert spaces 

c n ®c m = c nm 

e% ® e a — e{ a 
along with the corresponding identification of operator algebras: 

M n (C)®M ro (C) = M nm (C) 

e^j ® Cab — &ia,jb 

We denote by £, P, tt the objects associated in Definition 2.3 to h,k,h <g> k, with 
the choice of indices telling which is which. We have the following computation: 

(1) The ia-th row of h <g) k is given by (h ® k)i a = hi ® k a . 

(2) The corresponding magic basis is given by Ciajb = ® £,ab- 

(3) The corresponding magic unitary is given by Piajb = Pij ® Pab- 
Consider now the factorizations associated to h, k: 

A(S n ) -> A(h) - M n (C) 

A(5 m ) -» A(fc) -» M n (C) 
These are given by the following formulae: 

Xij > Uij > Pij 
•Eab * ^ab ¥ Pab 

We can form the tensor product of these maps, then we add to the picture the 
factorization associated to h®k, along with two vertical arrows: 

A{S n ) ® A(S m ) -» A(h)®A(k) M n (C)®M m (C) 
T I 

A{S nm ) -> A(/»®fc) M nm (C) 

Here the arrow on the right is the identification mentioned in the beginning of 
the proof. As for the arrow on the left, this is constructed as follows. Recall first 
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that A(S n m) is the universal algebra generated by the entries of a nm x nm magic 
unitary matrix Xj a j6- The matrix Xij <8) x a b being a magic unitary as well, we get a 
morphism of algebras mapping Xj a j6 — > <8) x a b, that we put at left. 
At level of generators, we have the following diagram: 



This diagram commutes, so the above diagram of algebras commutes as well. Thus 
the factorization associated to h®k factorizes through the algebra A(h)®A(k). Now 
since A(h ® k) is the minimal algebra producing such a factorization, this gives a 
morphism A(h) (g) A(k) — > A(h <g> k), mapping Uij <8> u a b — > Uj a j{,. 

This morphism is surjective, because the elements Ui a jb generate the algebra 
A{h®k). Now since A{h®k) makes a factorization of the representation associated 
toh(&k, this algebra produces as well factorizations of the representations associated 
to h, k. Thus our morphism is injective, and this gives the result. □ 

As a first application, we can start classification work for small Hadamard ma- 
trices. Recall from previous section that each matrix has a size n and a level I, 
according to the formula h € Hi(n) with I minimal, and that we decided to list 
matrices in terms of n, I, with n coming first. The situation is as follows: 

(1) For n = 1,2,3 there is only one matrix, namely the Fourier one. As shown 
by Theorem 2.1, this matrix produces the algebra C(Z n ). 

(2) For n = 4 the smallest possible level is I = 2. We have here the matrices 
M q , with parameter q = ±1. 

(3) The next possible level is I = 4. We have here the matrices M q with q = ±i, 
both equivalent to the Fourier matrix, which gives C(Z4). 

Thus at nl < 44, we have just two matrices left. But these can be investigated by 
using Theorem 3.1: 

Corollary 3.1. For q = ±1, the Hadamard matrix 



X{j X a b 

T 



'ia,jb 



Uij <g) U ab -> Pij® P a b 
I 

Uia,jb Piajb 



M q = 



(11 1 1\ 

1 q -1 -q 
1-11-1 

\1 -q -1 q J 



produces the algebra C(Z,2 x Z2). 

Proof. The 2x2 Fourier matrix is given by i^j = (— l) y , so we have: 




Consider the tensor square of F, given by (F (g) F)i a jb = F,i a Fjb- As for F, this 
is a ±1 matrix. Now since the —1 entry of F appears on the 11 position, the — 1 
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entries of F ® F appear on positions (ia,jb) satisfying ij = 11 or ab = 11, with the 
position (11, 11) excluded. There are six such positions: 

(11, 12), (11, 21), (12, 11), (12, 12), (21, 11), (21, 21) 

Now by using the lexicographic identification {1,2,3,4} = {11,12,21,22}, these 
six positions are 12,13,21,22,31,33. This gives the following formula: 

-1 -1 1\ 
-1 1 1 
1 -1 1 
1 1 l) 

By permuting rows and columns of this matrix we can get both matrices in the 
statement. The result follows from Theorem 2.1 and Theorem 3.1, by using the 
canonical identification C(%2) ® 0(1*2) = 0(1*2 x Z2). □ 

4. Magic squares 

In this section and in next one we investigate the case when A is commutative. 
In this situation, we must have A = C(G), for a certain subg roup G C S n . 

We already know that A is commutative for the Fourier matrix, where we have 
G = l n , and for the matrix M q with q = ±1, where we have G = I2 x 1*2- In fact, 
these two examples are the only ones that have been investigated so far. This makes 
our general problem quite unclear, because we have no counterexample. 

Let us go back however to proof of Theorem 2.1. A careful examination shows 
that the proof relies on the fact that P is a circulant matrix. The idea will be to 
generalize this proof, with "circulant" replaced by "magic" . 

Let us first recall the following well-known definition: 

Definition 4.1. A magic square is a square matrix a, having as entries the numbers 
1, . . . , n, such that each row and each column of a is a permutation of 1, ... ,n. 

We assume that all magic squares are normalized, in the sense that the first row 
is (1, . . . , n), and the diagonal is (1, . . . , 1). 

In this definition, we use the fact that it is always possible to normalize a magic 
square: a permutation of the columns makes the first row (1, . . . , n), then a permu- 
tation of the rows makes the diagonal (1, . . . , 1). 

As a first example, we have the n x n circulant matrix o~(i,j) = j — i, with 
j — i G {1, . . . , n} taken modulo n. Here is another example: 

2 3 4\ 
1 4 3 
4 1 2 

3 2 1/ 

The magic squares produce magic unitaries, in the following way: 

Proposition 4.1. If E = (E\, . . . , E n ) is a partition of unity with rank one projec- 
tions and a = cr(i,j) is a magic square, then (E a )ij = E a uj\ is a magic unitary. 



/I 



F®F 



V 1 



a 



(I 

2 
3 

V 4 
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In other words, E a is obtained by putting E superscripts to elements of a. For 
instance the above 4x4 matrix gives: 



E a = 



tE x 


E2 


E3 


Ea\ 


E2 


Ei 


E A 


E3 


£3 


E A 


Ei 


E2 


\E 4 


E 3 


E 2 


Ei) 



In the above statement, the condition that a is normalized is not necessary; nor 
the fact that the projections Ei are of rank one. We make these assumptions for 
some technical reasons, to become clear later on. 

We denote by 01, . . . ,o~ n the rows of a magic square a, and we regard them as 
permutations of {1, . . . , n}. For instance, for the above 4x4 matrix we get: 



o x = 


(1 


2 


3 


4) 


0-2 = 


(2 


1 


4 


3) 


0-3 = 


(3 


4 


1 


2) 


(74 = 


(4 


3 


2 


1) 



With these notations, we have the following result about commutativity, which 
deals with a slightly more general situation than that of Hadamard matrices: 

Theorem 4.1. For a representation ir : A(S n ) — > M n (C) having the property that 
Pij = ir(xij) are rank one projections, the following are equivalent: 

(1) A n is commutative. 

(2) A n = C(G), for a certain subgroup G C S n . 

(3) P = E a , for a certain partition of unity E and magic matrix a. 

Moreover, in this situation G is the group generated by the rows of a. 

Proof. Assume that A = A w is commutative. The algebra Im(7r) being a quotient 
of A, it is commutative as well. Thus the projections P^ mutually commute. 

On the other hand, two rank one projections commute if and only if their images 
are equal, or orthogonal. Together with the magic unitarity of P, this shows that 
each row of P is a permutation of the first row. 

So, consider the first row of P, regarded as a partition of unity: 

E = (Ei, . . . , E n ) 

By the above remark, the i-th row of P must be of the following form: 

E^ = (E ai (i),. ■ • ,E a .( n )) 

Here o~i <G S n is a certain permutation. Now the formula a(i,j) = ai(j) defines a 
matrix a, which is magic. With the above E and this matrix a, we have P = E a . 

In other words, we have (1) =^ (3). Also, the last assertion implies (2), which in 
turn implies (1). So, it remains to prove that (3) implies the last assertion. 
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Consider the following diag ram, where X — {<7i, . . . , cr n } is the subset of S n 
formed by rows of a: 

A(S n ) - M n (C) 



I T 

C(S n ) — C{X) 
We define the connecting maps in the following way: 

(1) The map A(S n ) — > M n (C) is the representation 7r. 

(2) The map — > C{S n ) is the canonical one, given by Xij — > Xij- 

(3) The map C(S n ) — > C(X) is the transpose of X C S' n . 

(4) The map C(X) -» M n (C) is given by <5 CTfc -» £ fc . 

We compute the image of Xij D Y the third map: 

Im(xij) = Im(x{<7 | cr(j') = 0) 

= I CTfe(j) = 

Thus at level of generators, we have the following diagram: 

x ij ^o-(i,j) 

I T 



This diagram commutes, so the above diagram of algebras commutes as well. Now 
since C(S n ) is a Hopf algebra, the Hopf algebra A w we are looking for must be a 
quotient of it. In other words, we must have A = C(G), where G is a subgroup of 
S n . On the other hand, from minimality of A n we get that this algebra must be the 
minimal one containing C(X). Thus G is the group generated by X. □ 

As an illustration, we get new proofs for Theorem 2.1 and Corollary 3.1: 

(1) For the Fourier matrix we have P = E a , where a is the magic square given 
by cr(i,j) = j — i, and E is the first row of P. We have Oi = c l , where c is 
the cycle c(i) = i — 1, so we get the group {1, c, . . . , c n_1 } = Z n . 

(2) For the matrix M q with q = ±1 we have P = E a , where a is the 4x4 matrix 
in the beginning of this section, and E is the first row of P. We get the 
group formed by rows of a, namely {cri,a 2 , cr 3 , 04} ~ Z 2 x Z 2 . 

In the general case, the construction a — > n — > A w is not understood. Some 
comments in this sense are presented in the end of next section. 

5. The commutative case 

We have now all ingredients for answering the question raised in the beginning of 
previous section. The result is as follows: 
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Theorem 5.1. For an Hadamard matrix h, the following are equivalent: 

(1) A is commutative. 

(2) h is a tensor product of Fourier matrices. 

Proof. The statement is of course up to equivalence of Hadamard matrices. We use 
this fact several times in the proof, without special mention: that is, we agree to 
allow suitable permutations of rows and columns of h, as well as multiplication of 
rows and columns by scalars of modulus one. 

(2) =4* (1): this follows from Theorem 2.1 and Theorem 3.1. 

(1) => (2): assume that A is commutative. Theorem 4.1 tells us that we have 
P = E a , for a certain magic square a. Here E is a certain partition of unity, but 
since a is normalized, E must be the first row of P. 

We have Pij = P(hj/hi), so the condition P = E a becomes: 




By using Proposition 2.2. we can assume h\ = 1. This gives: 

Now two rank one projections being equal if and only if the corresponding vectors 
are proportional, we must have complex scalars Ay such that: 

h - \ h 

h — A ij 'Mm) 

Now remember that the row vectors hi have entries of modulus one, so they are 
elements of the group T n , where T is the unit circle. We denote by Hi € T n /T their 
images modulo constant vectors in T n . The above condition becomes: 

This relation shows that X = {Hi,...,H n } is stable by quotients. Since X 
contains as well the neutral element Hi, it must be a subgroup of T n /T. 

Now since X is an abelian group, it must be a product of cyclic groups, and we 
can proceed as follows: 

(1) Assume first that we have X ~ Z n . We can permute rows of h, as to have 
Hi = H l ~ l for any i, for a certain element H G T n /T. This gives hi = Aj/9* -1 for 
any i, for a certain element p G T n , and certain scalars Aj. Now since the vectors hi 
are orthogonal, their multiples p t ~ 1 must be orthogonal as well. 

Consider now the vector p. Its image in T n /T has order n, so by multiplying h 
by a suitable scalar, we may assume that we have p n = 1. In other words, we have 
p = (wi), where each Wi is a n-th root of unity. 

Now from the relation p 1 ^ 1 _L fp^ 1 for 1 < i < j < n we get that the sum of fe-th 
powers of the numbers Wi vanishes, for 1 < k < n. This shows that the product of 
degree one polynomials X — W{ is the polynomial X n — 1, so the set of numbers Wi 
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is the set of n-th roots of unity. Thus by permuting rows of h we can assume that 
we have Wi = w i , with w = e 2m / n . 

Now the relation hi = \ip t ~ 1 tells us that h is obtained from the Fourier matrix 
by multiplying rows by scalars. Thus h is equivalent to the Fourier matrix. 

(2) Assume now that we have X ~ Y x Z , for certain groups Y, Z. By replacing n 
by nm and single indices by double indices, we can assume that we have H{ a = KiL a , 
where Y = {K%, . . . , K n } and Z = {L\, . . . , L m } are subgroups of T nm /T. 

We take now arbitrary lifts k\, . . . , k n and l±, . . . , l m for the elements of Y, Z. The 
formula Hi a = K{L a becomes hi a = \i a kil a , for certain scalars Aj a . Now these vectors 
being orthogonal, by keeping a fixed we get that the vectors k{ are orthogonal, so k 
is an Hadamard matrix. Also, with i fixed we get that I is an Hadamard matrix. 

On the other hand, we can get rid of the scalars Ajj by multiplying rows of h by 
their inverses. Thus we have hi a = kil a , where k, I are Hadamard matrices. 

With notations from previous section, this gives h = k I. 

(3) We can conclude now by induction. Assume that the statement is proved for 
n < N, and let h be a Hadamard matrix of order N. In the case X ~ Z^r we can 
apply (1) and we are done. If not, we have X ~ Y X Z as in (2), so we get h = k®l. 
Now Theorem 3.1 gives A(h) = A(k) ®A(l), so both A(k),A(l) are commutative. 
Thus we can apply the induction assumption, and we are done. □ 

The permutation groups coming from arbitrary magic squares can be non-abelian. 
Consider for instance the following matrix: 



A 


2 


3 


4 




3 


1 


2 


5 


4 


4 


5 


1 


3 


2 


2 


4 


5 


1 


3 


V 5 


3 


4 


2 


l J 



The rows of this matrix generate the group S$, which cannot be obtained by using 
5x5 Hadamard matrices, because of jHj- In general, we don't know what are the 
groups which can appear from the magic square construction. 

The other question is about what happens to Theorem 4.1 when the rank one 
assumption is removed. Once again, we don't know the answer. 

6. Small matrices 

We recall from [H] that for n = 1, 2, 3, 5 the only Hadamard matrix is the Fourier 
one, and that at n = 4 we have only the matrices M q , with \q\ = 1: 

1 1 1\ 

q -1 -q 
-1 1 -1 
-q -1 q J 

By results in previous sections, for parameters q satisfying q 4 = 1, this matrix M q 
produces a commutative algebra: 



(I 
1 
1 

V 1 
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(1) For q = ±1 we get C{G), with G = Z 2 x Z 2 . 

(2) For g = ±?we get C(G), with G = Z 4 . 

Now remember that for a finite abelian group G, the algebra of complex functions 
C(G) is canonically isomorphic to the group algebra C*(G). The isomorphism is 
given by the discrete Fourier transform in the case G = Z n , and by a product of 
such transforms, or just by Pontrjagin duality, in the general case. 

We can therefore reformulate the above result in the following way: 

(1) For g = ilwe get C*(G), with G = Z 2 x Z 2 . 

(2) For q = ±iwe get C*(G), with G = Z 4 . 

In order to be generalized, this result has to reformulated one more time. We 
have the following equivalent statement, in terms of crossed products: 

(1) For g = ilwe get C*(G), with G = Z 2 x Z 2 . 

(2) For g = iiwe get C*(G), with G = Zi x Z 4 . 

Recall now that for a discrete group G, meaning a possibly infinite group, without 
topology on it, the group algebra G*(G) is a G*-algebra, obtained from the usual 
group algebra <C[G] by using a standard completion procedure. 

With these preliminaries and notations, we have the following result: 

Theorem 6.1. Let n be the order of q 2 , and for n < oo write n = 2 s m, with m odd. 
The matrix M q produces the algebra C*{G), where G is as follows: 

(1) For s = we have G = Z 2n x Z 2 . 

(2) For s = 1 we have G = Z n / 2 x Z 4 . 

(3) For s > 2 we have G = Z n x Z 4 . 

(4) For n = oo we have G = Z x Z 2 . 

Proof. Consider the second and the third row of M q : 



P 

£ 



= (1 q -1 -q) 
= (1-11 -1) 



We regard these vectors as elements of the algebra C 4 . The product of these 
vectors being the last row of M q , we have the following formula: 



M a 



P 

e 



By taking quotients of these vectors we get the corresponding magic basis: 



/ 1 


P e 


ep\ 


p- 1 


1 ep- 1 


e 


e 


ep 1 


P 


{ep- 1 


e P- 1 
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Each entry of this matrix is of the form p k or ep k , with k € Z. So, consider the 
vectors p\ = p k and p k __ = ep k . These are given by the following global formula: 

4 = (1 ±<? fc (-l) fc ±(-g) fc ) 

Now the orthogonal projection onto a vector (xj) being the matrix (xjXj), we have 
the following formula for the corresponding projection: 



/ 1 



P{p k ±) 



±q K 



±q~ 
1 



±(-q)- k \ 
(-l) k 
±q- k 

1 / 



(-l) k 
±(-q) k 
±(-q)- k 1 
\±(-q) k (~l) k ±q k 

The idea is to express this matrix as a suitable linear combination of elements of 
C/(4). We write n = 2 s m as in the statement, and we proceed as follows: 

Case s = . Here n is odd, and we have q 2n = 1, so q n = ±1. We can replace if 
necessary q by —q, as to have c/ n = 1. 

Consider the canonical subgroup Z2 x Z2 C f7(4), consisting of the identity matrix 
1, and of the following three matrices: 





f° 


1 














1 


o\ 












1\ 




1 






















1 










1 





a = 











1 


1 











7 = 





1 















1 


°) 






1 





V 




K 1 








0/ 



Consider also the following matrix, which is in J7(4) as well: 

/0 g \ 





g - u J 
We have the following formulae, valid for k odd: 





q 


























\0 q- K / 
We also have the following formulae, valid for k even: 



\ 

,-fc 







(7 



3 













\ 



,-fc 



^- ft y 

: /3 and ct 2 ™ 






\o 

In particular we have a n = (5 and o"^™ = 1. 
For k odd the number n + k is even, and we have: 

P(p k ± ) = 1 ± aa n+k - a n T aa 
= (l-a n )(lTaa k ) 



[0 








q- k \ 








q k 





( 










\q k 








/ 


even: 








/0 c 







\ 






















q -k 


\o 





q k 


/ 
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For k even the number n + k is odd, and we have: 

P{p k ± ) = 1 ± aa k + a n T aa n+k 
= (l+a n )(l±aa k ) 

This gives the following global formula: 

P(4) = (1 + (-l)*<x n )(l ± (-l) k aa k ) 

It is routine to check that a,a £ U (4) generate the dihedral group G = 1,2n x ^2- 
Consider now the following element of the abstract group algebra C*{G): 

P k = (l + (-l) k a n )(l ± (-l) fc QCJ fc ) 

From <r 2n = 1 we get that <r n is self-adjoint, and of square 1. From a 2 = 1 and 
(aa) 2 = 1 we get aao = a, hence a k aa k = a for any It follows that aa k is also 
self-adjoint, and of square 1. This shows that in the above formula, both elements 
in brackets are projections. Moreover, these two projections commute. It follows 
that each element P± is a projection. 

Now remember that the magic unitary associated to it has as entries elements of 
the form P(p±). By making the replacement P(p±) — > P± we get a square matrix 
over C*(G), all whose entries are projections. The sums on rows and columns being 
1, this is a magic unitary, and we get a factorization of ir through C*(G). 

Case n = go . What changes in the above proof is that a, a £ U(4) generate now 
the infinite dihedral group G = Z x Z2. The second part of the proof applies as well, 
with this modification, and gives the result. 

Case s = l . Here n/2 is odd. We have q n = —1, so q n l 2 = ±«, which gives 
{—iq) n l 2 = ±1. We replace if necessary q by — q, as to have (—iq) n / 2 = 1. 

Consider the canonical subgroup Z4 C U(4), consisting of the identity matrix 1, 
and of the following three matrices: 



5 = 





1 





o\ 









1 






(0 








1\ 








1 





5 2 = 











1 




1 




















1 




1 
















1 








V 








V 




^0 


1 









Vo 





1 


0/ 



Consider also the following matrix, which is in U (4) as well: 

f-q \ 
q- 1 
T ~ -q 

\ <ry 

Since n/2 is odd and {—iq) n / 2 = 1, we have (ir)™/ 2 = 1. 
We have the following formulae: 



5r k = 



( 


q k 





° ^ 








(-q) k 














q -k 


\(-q) k 








/ 
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V 







(-q) k 





o / 



This gives the following formula for the above projection: 



P(p k ±) 



1 ± 5r k + (-l) k 5 2 ± (-l) k 5 :i T 



= (1 + (-l) k 5 2 )(l±5T k ) 

Consider now the matrix v = it. We have zW 2 = 1, and it is routine to check 
that 5, v G U(4) generate the group G = Z n / 2 x Z4. Consider now the following 
element of the abstract group algebra C*(G), where r 



-iv: 



Pi 



± ^ (l + (-l) k 5 2 )(l±5r k ) 

From <5 4 = 1 and from t5t = —5 we get that each P± is a projection. 

Now remember that the magic unitary associated to tt has as entries elements of 
the form P(p±). By making the replacement P(p±) — > P± we get a square matrix 
over C*(G), all whose entries are projections. The sums on rows and columns being 
1, this is a magic unitary, and we get a factorization of it through C*(G). 



Case s > 2 . Here n/2 is even, and we have q n = — 1, so q 



n/2 



We use the matrices 5, r. The above formulae for 5r k and <5 3 r fc hold again, and 
lead to the above formula for P(p±). What changes is the n/2-th power of t: 





f±i 





\ 




T n/2 _ 





=R 










±i 









V° 





TiJ 




we have r n = - 


-1, so 


{ W T) n = 


1, where w 



This shows that we don't have a factorization like in the case where n/2 is odd, 
so we must use the group generated by 5, wt £ U(4), which is G = Z n x Z4. □ 

The above result has several interpretations. First, it provides examples of sub- 
groups of the 4-th quantum permutation group, also called Pauli quantum group 
|BC2j . In other words, the whole thing should be regarded as being part of a gen- 
eral McKay correspondence for the Pauli quantum group, not available yet. 

Another interpretation is in terms of subfactors: the above result can be regarded 
as being a particular case of the ADE classification of index 4 subfactors. Unlike 
in the quantum group case, the classification is available here in full generality. See 
|EK| . The principal graphs corresponding to the matrices M q are those of type 

A4 , Dn \ Doo- This follows for instance by combining Theorem 6.1 with general 
results in |Ba2j . and with the well-known fact that cyclic groups correspond to type 
A graphs, and dihedral groups correspond to type D graphs. 

As a last remark, the above result provides the first example of a deformation 
situation for quantum permutation groups. This is of course just an example, but 
the general fact that it suggests would be use of the unit circle as parameter space, 
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for more general deformation situations. We should mention that this idea, while 
being fundamental in most theories emerging from Drinfeld's original work [D], is 
quite new in the compact quantum group area, where the deformation parameter 
traditionally belongs to the real line. It is our hope that further developments of the 
subject will be of use in connection with several problems, regarding both Hadamard 
matrices and compact quantum groups. 
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